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Abstract. 

We propose a jet model for the low/hard state of galactic black- hole X-ray sources which explains a) the X-ray 
spectra, b) the time-lag spectra, c) the increase of the amplitude of variability (QPO and high frequency) with 
■ increasing photon energy, and d) the narrowing of the autocorrelation function with increasing photon energy. 

Q\ ' The model (in its simplest form) assumes that i) there is a uniform magnetic field along the axis of the jet, ii) 

, the electron density in the jet is inversely proportional to distance and iii) the jet is "hotter" near its center than 

at its periphery. We have performed Monte Carlo simulations of Compton upscattering of soft photons from the 
accretion disk and have found power-law high-energy spectra with photon-number index in the range 1.5 - 2, 
power-law time lags versus Fourier frequency with index ~ 0.8, and an increase of the rms amplitude of variability 
and a narrowing of the autocorrelation function with photon energy as they have been observed in Cygnus X-l. 
Similar energy spectra and time variability have been observed in neutron-star systems in the island state. It is 
f*^ . therefore quite likely that a similar model holds for these sources as well. 
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I 1. Introduction with the non-thermal tail dominant and the rms variabil- 

V~J ■ ity of the source high (20% - 50%), then the source is in 

[Xj Galactic black-hole binaries have been extensively ob- thc i ow /hard state (see e.g. van der Klis 1995). 
^ \ served for more than thirty years. Their spectral and tim- 

■ - - ■ ing behavior exhibits a great variety of properties, yet in The location of thc energetic electrons that 

most cases they can be classified to be in one of few dif- Comptonize a soft photon input (and result in thc 

ferent states. The X-ray spectra of black-hole candidates power-law emission) is model dependent. One possibility 

have two distinct components: a multi-color black-body is that hot > optically thin plasma exists in a corona 

component (Mitsuda et al. 1984, but see Merloni et al. (Galeev et al. 1979; Haardt & Maraschi 1993; Haardt et 

2000) , widely accepted to come from the accretion disk aL 1994 ; Stern et aL 1995 ! Poutanen & Fabian 1999) that 
(Shakura & Sunyaev 1973), and a non-thermal hard X- lies above and below the thin disk - A second possibility 
ray tail that is usually modeled to come from unsaturated is a two-temperature flow (Shapiro et al. 1976; Ichimaru 
inverse Compton scattering of soft photons with energetic 1977 ! Rees et aL 1982 ! Narayan & Yi 1994; Abramowicz 
electrons (see e.g. Sunyaev & Titarchuk 1980; Poutanen et aL 1995 ! Esin et aL 1998 )- 

2001) . The relative strength of the thermal component in xhere is growing evidence that all X-ray binaries har- 
comparison with thc hard tail, the slope of the tail, and boring a black hole display radio emission when they are in 
the variability are used to identify the spectral state of thc i ow /hard X-ray state (Hjellming & Han 1995; Fender 
the source at a given epoch. When the luminosity of the et al 1997; Mirabel et al. 1998; Fender 2001; Stirling et 
source is low (less than a few percent of the Eddington lu- aL 2 001; Gallo et al. 2003). Inverse Compton scattering by 
minosity) and the thermal component is weak or absent, re l a tivistic electrons in a jet has been proposed as a mech- 

anism for the production of X-rays and 7-rays in Active 

Send offprint requests to: dgiannios@physics.uoc.gr Galactic Nuclei (Begelman & Sikora 1987; Bednarek et 
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al. 1996; Harris & Krawczynski 2002) and X-ray bina- 
ries (Band & Grindlay 1986; Levinson & Blandford 1996; 
Georganopoulos et al. 2002; Romero et al. 2002). The con- 
tribution of the synchrotron emission from the jet to the 
hard X-rays might also be significant (Markoff et al. 2001; 
Corbel & Fender 2002). 

Spectral information alone, however, can not reveal 
the accretion geometry. This is so because rather different 
models can result in very similar spectra. So as to make 
further progress in our understanding of how accretion 
works, studies of the timing properties of these sources 
are also essential. 

The fact that phase lags between two different X-ray 
bands are approximately independent of Fourier frequency 
(Miyamoto et al. 1988) is impossible to explain with mod- 
els where the Comptonization is taking place in a hot, 
uniform, optically thin plasma. This led to the idea of 
a large hot Comptonizing medium with density inversely 
proportional to radius (Hua et al. 1997; Kazanas et al. 
1997; Kazanas & Hua 1999; Hua et al. 1999). Still this 
model failed to explain the fact that the autocorrelation 
function of the light curve of Cygnus X-l in the low/hard 
state becomes narrower as the photon energy increases 
(Maccarone et al. 2000). 

Recently, Reig et al. (2003) proposed that inverse 
Compton scattering of soft disk photons by energetic elec- 
trons in a jet can explain both the X-ray emission of black- 
hole candidates in the low/hard state and the decreasing 
time lags with Fourier frequency. The density of the elec- 
trons in this (quite simple) jet model is assumed to drop 
inversely proportional of the distance z from the black 
hole. In such a case, and for an effective optical depth to 
electron scattering of order unity along the axis of the jet, 
there is equal probability for scattering of the photons per 
decade of distance z. The smallest period of variability 
that can pass through the scattering medium is of the or- 
der of the timelag caused by scattering, which is of the 
order of z/c. Faster variability gets smoothed out. This 
means that the maximum frequency of variability that 
can pass through the jet is inversely proportional to the 
timelag. Thus, timelags roughly inversely proportional to 
Fourier frequency can be produced. Applying this model 
to the case of Cygnus X-l it was possible to reproduce 
both the energy and the time-lag spectra. 

In this work we make one simple modification to the 
model described by Reig et al. (2003). Namely, we assume 
that the electrons are more energetic at small polar dis- 
tances p (in the core of the jet) than at larger p (jet's pe- 
riphery). With this modification we are able to explain two 
additional observational facts exhibited by Cygnus X-l. 
One is the narrowing of the autocorrelation function with 
increasing photon energy (Maccarone et al. 2000) and the 
other is the hardening of the high-frequency power spec- 
trum with increasing photon energy (Nowak et al. 1999). 
Preliminary results of this work appeared in Kylafis et al. 
(2004). A similar approach was followed by Lehr et al. 
(2000) to explain the energy dependence of the fractional 
amplitude of the high-frequency QPOs in GRS 1915+105. 



2. The model 

2.1. Characteristics of the jet 

Our jet model is quite similar to the one described in 
Reig et al. (2003). More specifically, we assume that the 
jet is accelerated close to its launching area and main- 
tains a constant (mildly relativistic) velocity thereafter. 
Furthermore, and for computational simplicity, we assume 
that a uniform magnetic field exists along the axis of the 
jet (taken to be the z axis). The electrons' velocity has 
two components one parallel (v\\) and one perpendicular 
(i>j_) with respect to the magnetic field. 

The new ingredient of the model is that the perpen- 
dicular component of the velocity is not constant, but it 
drops linearly with polar distance p, with its value in the 
outer edge of the jet being a fraction (typically half) of 
that along the z axis. The result is a jet whose core is 
"hotter" than its periphery. This, rather reasonable, as- 
sumption will have important implications on the shapes 
of the autocorrelations and the power-density spectra of 
the light curves in different energy bands. 

The electron-density profile in the jet is taken to be of 
the form 

n e (z) = n (z /z) p 1 (1) 

where z is the vertical distance from the black hole and p is 
a parameter close to 1 (see Sect. 3 below). The parameters 
no and zq are the electron density and the height at the 
base of the jet respectively. If H is the extent of the jet 
along the z axis, then for p = 1 the Thomson optical depth 
along the axis of the jet is 

T|| = no(T T z MH/zo)- (2) 

Let nR 2 (z) be the cross-sectional area of the jet as a 
function of height z. Then from the continuity equation 
M = irR 2 (z)m p n e (z)v\\ , we obtain the radial extent R of 
the jet as a function of height z 

R(z) = (i? 2 z/^o) 1/2 . (3) 

Here R = {M /-Km p noV\\ ) x l 2 is the radius of the jet at its 
base, M is the mass-ejection rate and m p is the proton 
mass. 

2.2. The Monte Carlo simulations 

For our Monte Carlo simulations we follow Cashwell & 
Everett (1959) and Pozdnyakov et al. (1983). A similar 
to ours, albeit independent, code was used in the work 
of Reig et al. (2001, 2003). The two codes are in excel- 
lent agreement between them. The procedure we follow is 
described briefly below. 

Photons from the accretion disk, in the form of a black- 
body distribution of characteristic temperature Tbb, are 
injected at the base of the jet with an upward isotropic 
distribution. It is not clear how the blackbody tempera- 
ture of the disk varies with polar distance p, since we do 
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not have a pure accretion disk, but we also have ejection 
of matter from the disk in the form of a jet. Thus, for 
simplicity, we consider a constant Tbb and constant lumi- 
nosity per unit polar distance p. These assumptions are 
not crucial for our conclusions. We also consider disk pho- 
tons that are emitted at poloidal distances p < Rq. This is 
justified by the fact that, for the typical values of Rq con- 
sidered here (~ 10 s cm, see Sect. 3), the disk flux outside 
this radius is small and furthermore, only a small fraction 
of it will interact with the jet. 

Each photon is given a weight equal to unity when it 
leaves the accretion disk and its time of flight is set equal 
to zero. The optical depth r(n) along the photon's direc- 
tion n is computed and a fraction e -r '") of the photon's 
weight escapes and is recorded. The rest of the weight of 
the photon gets scattered in the jet. If the effective op- 
tical depth in the jet is significant (i.e., >1), then a pro- 
gressively smaller and smaller weight of the photon expe- 
riences more and more scatterings. When the remaining 
weight in a photon becomes less than a small number (typ- 
ically 10 -8 ), we start with a new photon. 

The time of flight of a random walking photon (or more 
accurately of its remaining weight) gets updated at every 
scattering by adding the last distance traveled divided by 
the speed of light. For the escaping weight along a travel 
direction we add an extra time of flight outside the jet 
in order to bring in step all the photons (or fractions of 
them) that escape in a given direction from different points 
of the "surface" of the jet. The more a fractional pho- 
ton stays in the jet, the more energy it gains on average 
from the circular motion (i.e., v±) of the electrons. Such 
Comptonization can occur everywhere in the jet. Yet, a 
photon that random walks high up in the jet has a larger 
time of flight than a photon that random walks near the 
bottom of the jet. 

If the defining parameters of a photon (position, direc- 
tion, energy, weight, and time of flight) at each stage of its 
flight are computed, then we can determine the spectrum 
of the radiation emerging from the scattering medium and 
the time delay of each escaping fractional photon. The op- 
tical depth to electron scattering, the energy shift, and the 
new direction of the photons after scattering are computed 
using the corresponding relativistic expressions. 

2.3. The high-frequency power spectrum 

The power spectra of Cygnus X-l in the low/hard state 
can be well fitted by a doubly broken power low (Churazov 
et al. [2001]; see, however, Nowak [2000], Belloni et al. 
[2002] and Pottschmidt et al. [2003] for a different ap- 
proach and Revnivtsev et al. [2000] for evidence of a break 
above ~ 40 Hz). Other black- hole candidates show simi- 
lar power spectra when they are in the low/hard state 
(Miyamoto et al. 1992). The first break is typically ob- 
served at frequencies of the order of 0.1 Hz, where the 
slope changes from zero to ~ — 1 (when the power is ex- 
pressed in units of rms 2 /Hz) and the second break at 



about 1 Hz with the high Fourier-frequency slope being 
1.5. 

Nowak et al. (1999) showed, for the case of Cygnus X- 
1, that while the slopes of the first two power laws (slopes 
~ and ~ —1) do not depend on photon energy, the 
high-frequency power law exhibits a characteristic hard- 
ening with increasing photon energy. In this work we show 
that this hardening can be explained by a jet that has a 
core "hotter" than its periphery. Note, however, that other 
black-hole candidates show an energy independent power- 
density spectral shape (Lin et al. 1999). In the context of 
our model, this means that the jets in these sources have 
v±(p) w constant. 

In order to compute the power spectrum of the light 
curves in different energy bands, further assumptions con- 
cerning the source of the X-ray variability have to be 
made. We assume that the variability originates in the 
accretion disk, which is the source of the soft-photon in- 
put. This variability is modified as the disk photons travel 
through the jet and create the power-law energy spectrum. 

We think that it is reasonable to assume that the high- 
est frequency of variability, produced in a ring of disk ma- 
terial at radius p, is of the order of the Keplerian frequency 
at this radius, namely 

1*00 = 1/(2tt)(GM/p 3 ) 1 /2, ( 4 ) 

where M is the mass of the black hole. This assumption is 
justified by the small number of characteristic frequencies 
that describe the power spectra (Belloni et al. 2002). This 
then implies that a Fourier frequency v is associated with 
variability in the part of the accretion disk with radius 

p< (2ttv)- 2 / 3 (GM) 1 / 3 . (5) 

On the other hand, we have assumed that the core of 
the jet is "hotter" , resulting in a harder energy spectrum 
for input soft photons that come from the region close 
to the black hole. This is because these photons sample 
the core of the jet more than its periphery. Thus, higher 
Fourier frequencies are naturally associated with harder 
energy spectra, making the power spectra flatter in the 
harder energy bands. This is exactly what is observed in 
the case of Cygnus X-l (Nowak et al. 1999) and it will 
be quantitatively presented in Sect. 14.31 Furthermore, in 
the context of this model, the break of the power spec- 
trum that occurs at a frequency of a few Hz is naturally 
attributed to the Keplerian frequency at the outer edge 
p = Ro of the base of the jet. 

3. Parameter values 

Our jet model has a number of free parameters. The values 
of these parameters are chosen so that they are reasonable 
and they produce results similar to the ones observed. 

The parameters of our model are the following: The 
exponent p (Eq. [1]), the temperature of the soft pho- 
ton input, the extent H of the jet, the radius Rq and 
the height zq of the base of the jet, the Thomson optical 
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depth Tii along the axis of the jet, the parallel component 
of the velocity v\\ of the electrons, the perpendicular ve- 
locity component vx o at P — Q> an d its value at the outer 
edge (p = R) of the jet Ux,out- We assume that v± drops 
linearly with radius p, but any other smooth drop gives 
qualitatively similar results. We remark here that, due to 
the mildly relativistic flow in the jet, the emerging X-ray 
spectra from the jet depend on the angle between the line 
of sight and the magnetic field direction. 

The values of the parameters that reproduce quite well 
the observations of Cygnus X-l are here called reference 
values and they are: p = 1, fcTbb = 0.2 keV, H = 10 5 r„, 



Rq — 
0.55c 



100r„, z Q = 5r„, tii 



2vi 



= 2.5, uii = 0.8c, and v±,o = 
'-M.outj where r g = GM/c 2 — 1.5 x 10 6 cm stands 
for the gravitational radius of a 10 solar-mass black hole. 

In order to have good statistics in our Monte Carlo 
results, we combined all the escaping photons with direc- 
tional cosine with respect to the axis of the jet in the 
range 0.2 < cos6> < 0.6. Strictly speaking, this is incon- 
sistent with our assumption of a uniform magnetic field. 
If the magnetic field were uniform, then only one value 
of cos 8 would be appropriate for each source. As we will 
see (Sect. 4), a range of directional cosines is absolutely 
necessary in order to reproduce the observed phenomena. 
Therefore, a non-uniform (probably a spiraling) magnetic 
field permeates the jet. 

Besides the reference values, we have also explored 
what range of the parameter space is compatible with 
observations. Fortunately, this range is rather large. To 
quantify the last statement, we have found that, a varia- 
tion by a factor of two of parameters such as H, n , Rq, Zq 
and/or by 10% of the velocity components of the electrons 
or p, does not alter our conclusions. 

4. Results and discussion 

In what follows, we present results for the energy and time- 
lag spectra, the autocorrelation and cross correlation func- 
tions, and the power-density spectra. 




Fig. 1. Emergent photon- number spectrum multiplied by 
E 2 from our jet model for the reference values of the pa- 
rameters given in Sect. 3. The soft X-ray excess below 1 
keV is essentially the unscattered soft-photon input. When 
fitted with a power law with an exponential cutoff (see Eq. 
6), we find V = 1.7 and E cut = 280 keV. 




4.1. Energy and time-lag spectra 

The energy spectrum of the emerging photons for the ref- 
erence values of the parameters (previous section) is shown 
in Fig. ^ The spectrum above ~ 2 keV can be well fitted 
with a power law with an exponential cutoff 

dN/dE cx iT r e- E/Ecut , (6) 

where V — 1.7 and E cut = 280 keV. One can also see the 
soft X-ray spectrum (below ~ 2 keV), which is essentially 
the unscattered part of the soft-photon input. 

Keeping the rest of the parameters to their reference 
values and varying the Thomson optical depth along the 
axis of the jet T|| , we explore the role of this parameter 
on the emerging energy spectrum (Fig. 0). As expected, 
higher optical depths result in harder spectra (the pho- 
tons are scattered more and gain more energy from the 



Fig. 2. Emergent photon-number spectra multiplied by 
E 2 from our jet model for different optical depths tii . The 
dashed, dash-dotted, solid, dotted and dash-double-dotted 
curves correspond to tii = 1.5, 2, 2.5, 3.5, and 4 respec- 
tively. The rest of the parameters are kept to their refer- 
ence values. The power-law slope decreases with increasing 
optical depth from r ~ 2 (for t\\ = 1.5) to T ~ 1.45 (for 
tii = 4). The cutoff energy E cu ^ varies slowly in the re- 
gion 260 - 300 keV. For all the curves, the same number 
of input photons N — Ax 10 6 was used. 

energetic electrons of the jet). The photon number spec- 
tral index T decreases from 2 to 1.45, when tii varies in 
the range 1.5 — 4. 

We also made a similar study of the role of the radius 
Rq of the base of the jet. By increasing Rq, the spectrum 
hardens and vice- versa (Fig.^J). This is so because a larger 
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Fig. 3. Emergent photon-number spectra multiplied by 
E 2 from our jet model for different radii of the base of 
the jet. Here the dashed, dash-dotted, solid, dotted, and 
dash-double-dotted curves correspond to Ro/r g = 30, 50, 
100, 200, and 300 respectively. The rest of the parameters 
are kept to their reference values. The spectra harden as 
i?o increases, because the optical depth to electron scat- 
tering increases along the direction perpendicular to the 
axis of the jet. For all the curves, the same number of 
input photons N = 4 x 10 6 was used. 



R,q translates to a larger Thomson optical depth perpen- 
dicular to the jet axis (we kept the rest of the parameters 
at their reference values) and therefore more photon scat- 
terings. 

The cutoff energy E cut has a weak dependence on r\\ 
and Rq, but a rather strong dependence on v±. This is 
because v± determines the maximum energy gain of the 
soft photons. 

We then turn our attention to the timing proper- 
ties that can be extracted from our model and compare 
them with the observations. We start with the frequency- 
dependent time lags between two different energy bands. 
So as to simulate them, the time of flight of all escaping 
photons is recorded in 16384 time bins of duration 1/256 s 
each. The time of flight is computed by adding up the path 
lengths traveled by each photon and dividing the sum by 
the speed of light. Then we consider the light curves of two 
energy bands, 0.2 - 4 keV and 14 - 45 keV, so as to compare 
our results with observations directly (e.g. Nowak et al. 
1999). Following Vaughan & Nowak (1997), we compute 
the phase lags and through them the time lags between 
the two energy bands as a function of Fourier frequency. 

Figure 0] shows the time lags as a function of Fourier 
frequency for the fixed values of the parameters (Sect. 3). 
The time-lag spectrum is approximately represented by a 
power law with index /3 = 0.8 in the frequency range 
0.2 - 20 Hz, in good agreement with the observations of 
Cyg X-l and GX 339-4 (Nowak et al. 1999; Ford et al. 
1999). The low- frequency (~ 0.4 Hz) break of the time 



Fig. 4. Computed time lags as a function of Fourier fre- 
quency when the reference values of the parameters are 
used. Positive lags indicate that hard photons lag soft pho- 
tons. The results are logarithmically rebinned for clarity. 
The best-fit power law (in the 0.2 - 20 Hz frequency range) 
is i/~ - 8 . 



lags is determined by the range of directional cosines with 
respect to the magnetic field direction in the jet that we 
consider (Sect. 3). A larger range would produce a break 
at lower Fourier frequencies and vice versa. 

The time lags depend (among other things) on the 
dimensions of the jet. The scale of the jet is set by Rq 
and zq (see eq. 3). So as to quantify this dependence, we 
consider a "large" jet (Ro = 300r g , z = 20r g ) and a 
scaled down one (Rq = 50r g , z — 2r g ). In Figure [S] the 
time lags are presented for these two parameter sets and 
it is obvious that the larger jet results in larger lags by a 
factor of 2 for the same frequency. Notice, however, that 
the frequency dependence of the lags is very similar. 

Pottschmidt et al. (2000) studied the evolution of the 
lags for different epochs of Cygnus X-l and noticed that 
the lags have the tendency to increase in magnitude during 
state transitions. In the context of our model, this can be 
understood with the association of the state transitions 
with a more extended jet. 

4.2. Auto- and cross-correlation functions 

Maccarone et al. (2000) computed the autocorrelation 
functions (ACF) of the light curves of Cygnus X-l in differ- 
ent energy bands and showed that the width of the ACF 
decreases with photon energy. This imposes strong con- 
straints on the different models that attempt to describe 
both the spectral and the temporal properties of accreting 
black-hole binaries. 

With the simulated light curves of our jet model at 
hand (for the reference values of the parameters given 
in Sect. 3), it is straightforward to calculate the ACF at 
two different energy bands. The bands are kept to 0.2 
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Fig. 5. Computed time lags as a function of Fourier fre- 
quency for different dimensions of the jet. The crosses cor- 
respond to a large jet with Rq — 300r g and Zq = 20r g , 
while the open circles to a smaller jet with Rq = 50r ff 
and zq — 2r g . The lags of the bigger jet are a factor of 2 
larger than these of the smaller one at the same frequency. 
In either case, the best-fit power law (in the 0.2 - 20 Hz 
frequency range) is z/~ ' 8 . 



Fig. 6. The autocorrelation functions of the simulated 
light curves in two energy bands is shown. The thin solid 
curve corresponds to the soft band (0.2 - 4 keV), while 
the thick solid curve to the hard band (14 - 45 keV). 
Both curves are normalized to 1 for lag r = 0. Notice 
that the autocorrelation function becomes narrower in the 
higher energy band. We also plot (dashed curve) the cross- 
correlation function of the two light curves. 



- 4 keV (soft band) and 14 - 45 keV (hard band). The 
ACFs are normalized to unity at zero lag in both bands 
and are shown in Fig. We find that the width of the 
ACF decreases with photon energy in qualitative agree- 
ment with observations. At first sight, this is contrary to 
intuition. It can never happen in a spherical Comptonizing 
medium. However, in our model the high-energy photons 
are not escaping isotropically, i.e., with a uniform distribu- 
tion in cos 8. The more energetic photons come out pref- 
erentially at large values of cos 9 and with a lightcurve 
which is narrower than that of the less energetic photons. 
This gives rise to a narrowing of the ACF with increasing 
photon energy. In a realistic model with a non-uniform 
magnetic field, this could come about as follows: High- 
energy photons (escaping with large values of cos 9) come 
from relatively small regions of the jet and have a nar- 
rower lightcurve than the softer photons which come from 
the entire jet. 

In Fig. El we also plot the cross-correlation function of 
the two light curves. It should be stressed that the exact 
shape of the auto- and cross-correlation functions depends 
on the properties of the shots of the soft-photon input. 

4.3. Power-density spectra 

Following the method described in Sect. 12.31 we calculate 
the high-frequency part of the power spectrum in five dif- 
ferent photon energy bands (0.2 - 4 keV, 4-6.3 keV, 6.3 

- 9 keV, 9 - 14 keV, and 14 - 45 keV). 

In the context of our model, the second break in the 
power spectrum corresponds to the Keplerian frequency 



at p = Rq in the accretion disk, which for Rq — 100?* g = 
1.5 x 10 8 cm is z/br — 3 Hz. At higher frequencies, the 
rms soft-input variability is assumed to have a power-law 
frequency dependence 

( v \~ 1/3 ,«\ 
rmsi n = rms {T^J , (') 

so as to match the high-frequency slope of the power spec- 
trum for low-energy photons {E ~ 1 keV). Note that for 
rmsin oc i/ -1 / 3 we have rms 2 n /Hz oc z; -5 / 3 . Here rmso is 
the observed rms at 3 Hz. 

We remind the reader that, according to our assump- 
tions (Sect. 12. 3f) . the variability at a given Fourier fre- 
quency v is identified with the one that comes from the 
part of the accretion disk with r constrained by expres- 
sion (5). Thus, for a given Fourier frequency v, we vary the 
soft-photon input from the appropriate part of the accre- 
tion disk by an amount determined by Eq. (7) and record 
the resulting variability in the different energy bands. We 
find that, the higher the Fourier frequency the larger the 
fractional change in rms 2 /Hz between the soft band (0.2 
- 4 keV) and one of the higher bands (Fig.EJ. This can be 
understood as follows: The higher the Fourier frequency 
the closer to the core of the jet the source of variability is 
and thus the soft-input photons get upscattered in a "hot- 
ter" medium. By construction of our model, the curves of 
rms 2 /Hz for the different energy bands meet at ~ z^b r , 
which corresponds to r = Rq, i.e., the entire accretion 
disk considered in our model. 

As a consequence of the above, the high-frequency 
(y > z^, r ) slope of the power spectrum depends on the en- 
ergy band. A power-law fit to the power spectrum for the 
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Fig. 7. Power-density spectrum of Cyg X-l. The dashed 
line corresponds to the low-frequency part of the power 
spectrum, as typically observed. The part of the power 
spectrum above the second break is the result of our 
simulations and shows the spectrum's photon-energy de- 
pendence. The five curves correspond to different en- 
ergy bands as follows: b\ : (0.2 — 4), 62 : (4 — 6.3), 63 : 
(6.3-9), 64 : (9-14), 65 : (14-45) keV. The steepest slope 
(thinner curve) corresponds to b\, the second in steepness 
to &2 and so forth. 

soft energy band (0.2 - 4 keV) results in an index 1.7, es- 
sentially unchanged from the input value 5/3. As the band 
energy increases, the power spectrum flattens and for the 
hardest energy band considered (14 - 45 keV) the index 
becomes 1.4 (Fig. EJ , in agreement with the observations 
of Cyg X-l (Nowak et al. 1999). This variation of the index 
by ~ 0.3 is a consequence of the negative radial gradient 
of the perpendicular component of the electron velocity 
vj_ and depends on how much more energetic are the elec- 
trons in the core of the jet in comparison to those at its 
periphery. Here we have taken v± t o = 0.55c = 2ui^ ut- 
Had we taken v±^ = u_i_ l0 ut (as in Reig et al. 2003), no 
flattening of the power-density spectra with photon en- 
ergy would have been produced. This may be the case for 
some black- hole candidates (Lin et al. 1999). 

5. Conclusion 

The non-thermal X-ray emission in black-hole binaries is 
usually modeled to come from inverse Compton scattering 
of soft photons in a hot corona (e.g. Poutanen & Fabian 
1999) that lies above and below a thin disk (Shakura & 
Sunyaev 1973) or in a two-temperature flow (Shapiro et 
al. 1976). Here we explore an alternative picture where the 
inverse Compton scattering takes place in the jet whose 
presence, whenever the source is in the low/hard state, is 
well established (Fender 2001; Stirling et al. 2001; Gallo 
et al. 2003). 

In a previous work (Reig et al. 2003), it was shown 
that, for a range of the parameters of a very simple 



Fig. 8. Power-law slope of the high-frequency power- 
density spectra (Fig. [TJ) versus photon energy. The abo- 
lute value of the slope for the different energy bands is 
61 : 1.72, b 2 : 1.61, b 3 : 1.56, 64 : 1.52, 65 : 1.41. 

jet model, both the energy spectra and the frequency- 
dependent time lags of Cygnus X-l could be reproduced. 
In this work we have added a negative gradient along the 
polar direction r of the electron velocity component v±, 
resulting in a jet whose core is "hotter" than its periphery. 
This jet model can reproduce not only the energy spec- 
trum and the time lags but also the fact that the autocor- 
relation function is measured to be narrower (Macarrone 
et al. 2000) and the high-frequency power spectrum harder 
(Nowak et al. 1999) for the harder photon energy bands. 

In the context of this model the X-ray variability comes 
from the soft photon input that is further reprocessed in 
the jet. The high frequency break in the power spectrum 
(typically seen at ~1 Hz) is identified with the Keplerian 
frequency at the outer edge of the base of the jet. 

One important challenge to our model is to demon- 
strate that, with no change of the assumptions, one can ex- 
plain quantitatively the observed radio spectrum of black- 
hole sources in the low/hard state (Fender 2001) and the 
correlation between radio and X-ray luminosities (Gallo 
et al. 2003). Preliminary calculations indicate that this is 
possible. A detailed presentation is upcoming (Giannios, 
in preparation). 

Neutron-star systems in the island state exhibit energy 
spectra and time variability similar to those of black-hole 
systems. We therefore suggest that our model may be rel- 
evant for these neutron-star systems as well. 
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